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^ ■ Abstract 

^ I We consider supercritical accretion onto a compact object when the accretion 

^ ■ rate exceeds its value derived from the Eddington limit: MEdd = LEdd/c^- 

We use the scenario that was proposed by Shakura and Sunyaev (1973) for 
the supercritical regime, according to which matter flows out of the gas-disk 

^ ! surface under radiation pressure. An analytic solution of the accretion-disk 

^ ' structure is obtained for this scenario by using classical equations. An at- 

tempt is made to construct a model for nonconservative supercritical ac- 

^ ■ cretion with advection. To this end, we develop a numerical scheme for 

calculating the structure of an advective disk with mass outflow. The goal 
of solving the problem of a supercritical disk with mass loss is to explain 
accreting sources with outflowing envelopes, for example, SS 433 and galac- 
tic nuclei. The structure and observational manifestations of the envelope 
which forms around a supercritical disk are considered by taking into account 
bremsstrahlung absorption and Thomson scattering of photons. 



Introduction 

The accretion- disk theory is used in astrophysics to explain a wide range of X- 
ray sources, their energetics, spectra, and variabihty. By now a wide variety 
of models have been proposed to account for the observed characteristics of 
quasars. X-ray pulsars, and cataclysmic variables. 

The disk-accretion studies by Gorbatskii (1964), Lynden-Bell (1969), 
Shakura (1972), Pringle and Rees (1972), and Shakura and Sunyaev (1973) 
yielded the first mathematical theory of accretion disks. The general pattern 
of steady-state accretion in disks is as follows: the gravitational energy of 
the gas that rotates around the central star transforms into the thermal and 
kinetic energy of orbital motion through friction. For the effective emission 
of heat, the gas flows cool down and produce a disk structure. The disk lu- 
minosity is a function of the central object's mass, the rate of accretion onto 
it (i.e., the rate of mass inflow to the outer boundary of the accretion disk 
from the interstellar medium or from a neighboring star), and the specific 
mechanism of angular-momentum transfer in the disk. 

When the accretion rate is not too high, the total accretion-disk luminos- 
ity is directly proportional to the accretion rate. Under such quiet conditions, 
the global disk structure can be successfully described by the classical solu- 
tion of Shakura and Sunyaev (1973). 

However, it is well known from the scenario for the evolution of binaries 
(Tutukov and Yungelson 1973; Van den Heuvel 1994) that there are evolu- 
tionary stages of violent mass transfer between the component stars in binary 
systems. At M greater than the critical value 

Mcr = LEddMc^ ^ 2 X 10-^ - Mo/yr , (1) 

v 

where m = M/Mq is the stellar mass (in solar masses), and i? is the ac- 
cretion efficiency in a thin disk, the disk luminosity becomes so large that 
the radiation pressure begins to hamper the infall of matter to the center. 
This accretion regime is called supercritical or super-Eddington. In a spher- 
ically symmetric case, the luminosity reaches the Eddington limit, when the 
gravitational attraction by the central star equals the radiation pressure on 
a proton: 

LEdd '^T rrijjMG oo 
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LEdd ~ 1.25 X lO^^m erg/s , (2) 



where a = 6.65 x 10~^^cni^ is the cross section for Thomson scattering by 
free electrons. 

In reahty, the accretion rate clearly exceeds the critical value (0) in several 
cases. For example, in SS 433, matter from the supergiant flows onto the 
compact object (a neutron star or a black hole) at a rate of ~ 10~'^Mq yr~^ 
(Cherepashchuk 1981). High accretion rates may occur in galactic nuclei 
during intense infall of stellar matter onto supermassive black holes (see, 
e.g., Dopita 1997; Szuszkiewicz et al. 1996). 

Shakura and Sunyaev (1973) proposed a self-regulation mechanism for 
supercritical accretion. As soon as the luminosity of the disk exceeds the 
Eddington limit, matter flows out of its surface, causing the rate of accretion 
onto the central object to decrease. The amount of gravitational energy 
which can be liberated in the disk thus decreases, and the disk luminosity is 
limited by a value that does not greatly exceed the Eddington limit. 

The outflow from the disk surface takes place inside the radius at which 
the disk thickness is equal in order of magnitude to the disk radius. The 
disk structure outside can be described by the classical Shakura-Sunyaev 
solution. This radius Rg is called the spherization radius. As was established 
by Shakura and Sunyaev (1973), if the accretion rate in the disk decreases 
as 

M{r) = ^Mo , (3) 

then the total disk luminosity does not greatly exceed the Eddington limit: 
Ltot ~ LEdd X (1 + ln{Mo/Mcr)) ■ In this case, the total mass loss rate in 
the wind is ~ Mo(l — Mcr/Mg) , where Mo is the initial accretion rate at the 
outer disk radius. 

The hypothesis of mass outflow from the disk is confirmed, for example, 
by photometric and spectroscopic observations of SS 433, which suggest the 
presence of an envelope of matter flowing out of the disk (e.g., Murdin et 
al. 1980). In some galactic nuclei, observations reveal outflow rates reaching 
~ 1Mq/jt (Heckman et al. (1981). 

The authors who consider the winds from accretion disks, for example, in 
cataclysmic variables, propose the mechanism of wind generation in spectral 
lines. The rate of mass loss by such disks is ~ 10^"*^^^^^ Mq/jt at a flow 
rate of ~ 10~^ Mq/jt (see, e.g., Suleimanov 1995; Pereyra et al. 1997). By 
contrast, for supercritical mass outflow, a large optical depth in the disk and 
the wind gives rise to a blackbody spectrum and causes this mechanism to 



decrease in importance. Meier (1979) considered the wind generated by a 
supercritical disk at a certain radius close to the spherization radius (for the 
Shakura-Sunyaev disk model) and computed the vertical structure of the 
disk and the wind. Here, we consider the radial structure of an outflowing 
disk with height-averaged physical parameters. For SS 433, our outflow rate 
is comparable in order of magnitude to the value of Meier (1982): 3 4- 9 x 
10"^ MqI yr. 

Equations for disk accretion with mass loss 

The structure of accretion disks at subcritical accretion rates is described by 
standard equations of disk accretion (Shakura and Sunyaev 1973). In this 
case, the disk geometry can be roughly described by a thin disk, and ap- 
propriate simplifying assumptions, for example, replacement of the physical 
parameters in the disk by their height-averaged values, can be imposed on 
the equations. To describe a nonconservative disk — a disk with mass out- 
flow from its surface, we use the equations from the standard model, which 
we modified to fit our view of mass outflow from the disk. 

The equation of radial angular-momentum transfer along the disk — the 
ip component of the Euler equation integrated over the coordinate perpen- 
dicular to the disk plane, is 

l{M{rW] = -2.l{W,,r-] + ^ur^ , (4) 

Here, M{r) is the accretion rate, r is the radius, u is the Keplerian angular 
velocity, and Wr<p is the integrated, height-averaged tangential component 
of the stress tensor Wr^p in the accretion disk: 



Wr^ = 2HWr^ Wr^ = -Tjr -^ i (5) 



dr 



where H is the disk half-thickness, and rj is the dynamical viscosity. The 
change in angular momentum of the accreted gas is the sum of the angular- 
momentum transfer by viscosity (the first term on the right-hand side of 
(^) and the transfer of angular momentum carried away by the outflowing 
matter (second term). This equation differs from the standard equation of 



angular-moment um transfer (see Eq. (p!^), in which M = const, precisely 
by the presence of the second term on the right-hand side. 

The turbulent viscosity can be described by the dynamical coefficient, 
?7 ~ p < Vtk > , where Vt and It are the characteristic velocity and length of 
turbulent mixing, respectively. Using Prandtl's semiempirical relation 



we find from (^) that 



Vt = kr-r- 
dr 



Wrip — —pVt = —f^i PCt^ 



where rrit is the Mach number, mf = v^/o? ,, and a is the speed of sound. 
Shakura (1972) introduced a turbulence parameter a, which relates the tan- 
gential component of the stress tensor to the pressure in the disk P = 
Pgas + Prad (the sum of the gas and radiation pressures) : 

Wrip = —OiP . (6) 

The disk models that assume turbulence to be the source of viscosity and 
that use relation (|^) are called a-disks. 

The equation of energy balance: the mechanical energy due to viscosity 
produced in the disk at radius r is completely emitted at the same radius, 

g+(r) = Qrad{r) , (7) 

where Qrad{f) is the radiative fiux from one of the two accretion-disk surfaces, 

and Q~^{r) {Q~^{r) (erg/cm^s) is half the heat released in a disk annulus of 

width dr. 
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where k = a /rrip = 0.4 cm^/g is the specific Thomson scattering cross 
section (jUp = 1.67 x 10"^"^ g is the proton mass), Srad is the mean radia- 
tive energy density, and p is the mass density. Equation (|^) requires special 
stipulations in the case of a high accretion rate: it turns out that the heat 
exchange between a given portion of the disk and its adjacent portions must 
also be taken into account, as discussed below. Note that the "local" equa- 
tion of energy balance (|^) allows an analytic treatment of the problem with 
outfiow. 



Q+{r) = --uWr^ , Qrad = - ,^ ^jT ^rad , (8) 



The expression for Qrad in (||) is the solution of the transfer equation 
in the diffusion approximation, which is vahd for optically thick disks. The 
optical depth r of such a disk, whose geometrical thickness can be determined 
from the condition that the vertical component of the gravity is equal to the 
pressure (hydrostatic equilibrium), 

IdP P _ GMH 

is much larger than unity at supercritical accretion rates. Here, we disregard 
the general-relativity effects. Novikov and Thorne (1973) considered thin 
accretion disks around Schwarzschild black holes in terms of the general- 
relativity theory. 

The equation of state 

P = pa^ = ^I^ (10) 

o 

holds in the radiation-dominated regions of thin disks. In the supercriti- 
cal regime, the radiation pressure in the disk is much higher than the gas 
pressure. 

The spherization radius Rg at which mass outflow sets in can be ap- 
proximately determined from the condition that the luminosity reaches the 
Eddington limit (Lipunov 1992): 



^''^*'°- .,.. theuwehave R,.R.M^, (11) 



2 R, ^^'"'' ^^' ^"-"M., 



where Ro is the radius of the inner disk boundary. 

The matter flowing out of the disk surface produces an expanding en- 
velope. The problem of mass outflow from a supercritical accretion disk is 
similar, except for the geometry, to the problem of outflow of strong winds 
from WR stars and giants. The wind structure is described by the continu- 
ity, momentum, and energy equations. We do not set the goal of solving this 
complex problem and use an upper limit to determine the outflow rate. For 
the matter from the disk to be moved to infinity, the work dM (^^ ^^] 

must be done in a second. Clearly, an upper limit on dM can be placed by 
assuming that this power is equal to 2 Qrad'^'^rdr, the luminosity of the disk 
annulus from the two surfaces. We thus obtain the following estimate, which 



is the so-called "energy-wind" relation: 

/,)2t- rlM(r\ 

Qrad ■ (12) 



1 dM (vl - vl^^i) _ J'r dM{r) 



2 27rrrfr 2 Stt dr 

For the mass outflow from the disk, the accretion flux decreases with 
radius and is the sought-for function which is given by 

M{r) = AnrHpVr , (13) 

where Vr is the radial velocity in the disk. 



Conservative solution 

The solution to the set of equations (|) (without the second term on the 
right-hand side)- ([lO|) is a classical solution for a thin accretion disk without 
outflow, which is called the standard model. In this case, the accretion 
efficiency is 6' = 1/12. The equations can describe the disk structure in 
different regions, depending on the form of equation of state, i.e., on the 
relative contribution of the radiation pressure and the gas pressure. The 
solution to the equation of angular-momentum transfer, 

^^{Mu;r^} = -27r^^{Wr,r^} , (14) 



for M = const and the Keplerian law of motion lu = JGM/r^ is the following 
(Shakura and Sunyaev 1973) 



1^^- ^hM"^ 



Wr, = ^{i-\l^\+w:;{^) , (15) 



where Rin is the inner radius, and Wr^p{Rin) = W^^ . In the standard solution 
at the inner disk boundary _Rj„ = Ro, one assumes that W^^ = 0. 

The total disk luminosity in this solution is proportional to the accretion 
rate, 

Ltot = 7^—^- = LEdd -— , M <M,r . (16) 

/ rCo Mrr 



If the central star is a Schwarzschild black hole, then Ro = 3Rg, the 
radius of the marginally stable orbit; if it is a neutron star, then the inner 
disk boundary lies either on the neutron-star surface^] or at the Alfven radius 
in the case of a strong magnetic field near the neutron star. 

The disk thickness in the regions under consideration changes with radius 

only slightly, as we show below. The condition for radiation domination can 

be written as the equality of the forces of attraction and radiation pressure 

on the disk surface, 

mpGMH Qrad o- 

^^ ■ 17 

Using d^), (H), and (0), we obtain Qrad{,f) and, after simple algebraic trans- 
formations, derive 

M 3 f [r^\ M 3 , 

H = Rg^—\l-J—\^Rg^—{r»Ro). 18 

Solution in a nonconservative model 

The set of equations (^-(|13]) describes the disk structure within the spher- 
ization radius, where mass outflow proceeds. Its solution can be found ana- 
lytically. From Eqs. (§, (0), (|), and (|T2|), we derive M{r) and Wr^ : the 
accretion rate in the disk, starting from the spherization radius Rs at which 
outflow sets in, falls off to the center as 



R.\"'^ + l[i; 



r 



5/ 
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M{T)^MXr)( — ) ) i^n • ^<R.- (W) 



i + i(f) 
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This law has the asymptotic limit 



M{r)^Mo{r)^ Ro « r < R^ , 

Kg 

i.e., the derived solution matches Shakura-Sunyaev's solution (Eq. (|^)). The 
analytic expression for the height-averaged component of the stress tensor in 



"'^For a neutron star with a soft equation of state, the stellar radius may turn out to be 
smaller than the radius of the marginally stable orbit. 



the disk is 

^,,.^M!>W4) . (20) 

'' i + Kfe) 

The boundary conditions at the inner radius Ro are Wr^p = and Qrad = 0. 
The absence of stress at the inner disk edge is undoubtedly a simphfying 
assumption, which ignores the hydrodynamic pattern of the flow. 

The full solution is obtained by joining the solution with outflow and 
the solution without outflow (outside the spherization radius), which is the 
standard Shakura-Sunyaev solution for the appropriate boundary conditions 
at the spherization radius Rg (Eq. ([T5|) ). 

The spherization radius is a parameter of the problem. It can be numer- 
ically adjusted in such a way that the disk luminosity outside this radius is 
of the order of LEdd '■ 

Rs : L{r > R^, i?,) = LEdd • (21) 

Note that the refined spherization radius gives a different total rate of mass 
outflow from the disk. Bisnovatyi-Kogan and Blinnikov (1977) studied the 
dynamics of particle motion near a transcritical disk and found that the 
process of particle escape to infinity under radiation pressure in the gravita- 
tional field of a central object or disk spherization is established at a total 
disk luminosity in the range 0.6 to ILEdd- It was found that up to 80% of the 
matter could be carried away from the disk for O.GLEdd However, since Rg is 
currently rather uncertain, we use condition ( PT| ) in the numerical procedure. 

In Fig. |I|, the solid lines represent some of the disk characteristics as a 
function of radius for rh = M/Mcr = 1000; the change in accretion rate, 
which is described by Eq. (plQ]), is shown in the lower left panel; the derived 
spherization radius is Rg = 1.62 Ro Mo/Mcr (cf. Eq. (|TT|)). 

The upper left panel shows a plot of the relative half-thickness H/r in the 
disk against radius, which is constant (solid line) for r < Rg and is the same 
for any Mq > Mcr- This can be easily verified by substituting the accretion 
rate M{r), which decreases under the effect of the energy-wind mechanism 
(p!9D, into the expression for the disk thickness in the radiation-dominated 
region (pIS]). 



Advective disks with mass loss 

Formulation of the Problem and Results 

At supercritical accretion rates, the disks are always thick. The assumption 
that H/r << 1 is invalid for them. This can be illustrated by an analytical 
solution for a nonconservative disk: we see from Fig. |I| (upper left panel, 
solid line) that the relative half-thickness reaches ~ 0.6. 

In a thick disk, the emission from the surface is not an efficient cooling 
mechanism. In the characteristic time of inward matter motion along the 
radius, the heat cannot radiate away and is transported with the matter, 
i.e., the time it takes for a photon to reach the disk photosphere, tjjf = 
3Ht/c (here, r = Hna^ >> 1 is the optical depth for scattering across 
the disk, and n is the height-averaged number density of the disk matter), 
is much longer than the characteristic time of radial particle motion r /vr, 

where Vr ^ a{^\ ur is the radial velocity with which the matter moves 
over the disk (see Fig. 0). 

This transport of heat along the disk with matter is called advection. If 
advection is taken into account, the equation of energy balance (^ changes: 

Q = Qrad + Qadv , Qadv = pVrHT—— , (22) 

dr 

where T is the temperature, and 5* is the specific entropy. Using standard 
thermodynamic relations, the expression for advection can be rewritten as 



dU)+p d{^] 
Qad. = pVrH ^'^ ^^ ^ , (23) 

where e is the internal energy density of the gas, and P is the total gas 
pressure (averaged over H). Inside the disk, where the optical depth is 
much larger than unity, the photons and particles are in thermodynamic 
equilibrium, and e = AT'^, A = 7.58 x 10~^^erg/cm^^. Since the radiation 
pressure dominates over the gas pressure in supercritical disks, we have 

P ^ Prad = ^- . (24) 

Thus, the set of equations (|)-(|), (|), d), dH), (ll), and (H-dl) 



describes an outfiowing disk with advection. A numerical solution of the 

10 



disk structure can be found by using a developed computer code, which has 
two modifications: to compute the structure of an advective disk with and 
without outfiow. We tested the numerical scheme on the Shakura-Sunyaev 
standard-disk model. 

The results of our calculations for the accretion rate Mo/{12LEddl(?) = 
1000 are shown in Fig. |I| (long dashes). The short dashes represent the results 
of our calculations for a conservative advective disk: no mass loss by the disk, 
and the accretion rate is constant. 

As we see, in the solutions with advection, it dominates in the inner disk 
regions, i.e., only a small fraction of the heat is emitted: QradlQ'^ « 1- 
In this region, the disk cooling by the energy wind (Eq. (|12D) is inefficient. 
Nevertheless, it is clear that a considerable fraction of the initial accretion 
rate can be lost by the disk (long dashes). 

Our numerical calculations for various m in the range lO^-lO'' yielded 

i^W^Edd-O.e + O.rinm. (25) 

Figure ^ shows plots of the accretion rates in the disk for various values 
at the outer disk boundary. 

Observations of SS 433 reveal an expanding envelope around the accre- 
tion disk. The binary system SS 433 consists of an O or B star and a com- 
pact object (a black hole, M ~ lOM©, or a neutron star, M ~ IM0) (see, 
e.g.. Van den Heuvel 1981). Mass accretion onto the compact object on 
the thermal time scale of a normal star proceeds in the system at a rate of 
lO-^-lO-^Mo/yr (Cherepashchuk 1981; Van den Heuvel 1981). The mass 
loss rate in the envelope is lO~^-lO~^M0/yr (Shklovskii 1981; Dopita and 
Cherepashchuk 1981; Van den Heuvel 1981). 

Thus, the observed outfiow rate from the system, which is comparable to 
the accretion rate, is obtained in the model of a supercritical nonconservative 
disk. 

Envelope around a supercritical disk 

Luminosity and Effective Temperature 

The envelope that is formed from the matter fiowing out of a supercritical 
disk can be optically thick for the emission from the inner disk regions. Its X- 
ray emission is reprocessed into a softer emission. Let us roughly estimate the 
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effective temperature and emergent luminosity of tlie envelope as a function 
of the model parameters and the empirically derived quantities 

jylp„ , iio , Rn , . 

•^' = 14- ''==^' -"'=«;■ '''' 

The first two coefficients are determined when the disk structure is calculated 
numerically; Me„ is the mass loss rate in the envelope; M^ is the initial 
accretion rate at the outer disk boundary; Ro is the inner disk radius; Rg is 
the spherization radius which is derived from condition (|21|); and the third 
coefficient is equal to 3, if the compact accreting object is a black hole, and 
is set approximately equal to 3 in the case of a neutron star: 

Ax ~ 0.5 , A2 = 0.87 , ^3 = 3 . 

We assume that the mass outflow proceeds at a constant rate 

Men = AirRlmprisV^iRs) , (27) 

where Ug is the number density in the envelope at radius Rg, and foo(-Rs) = 



2GM/Rs = 0.62 c m~^/^ is the outflow rate. 
The emergent spectrum differs in shape from a blackbody spectrum, be- 
cause the scattering of emergent photons by free electrons plays a major role 
in the opacity (Shakura and Sunyaev 1973; Sunyaev and Shakura 1974). The 
visible envelope radius is R^, where the photon is scattered for the last time. 
The radius R„ can be derived from the relation 

/(J n R 
a^ndr = -^ = 1 , (28) 

R 

T 

by taking into account the fact that the particle number density in the en- 
velope at a constant outflow rate in it falls off as the inverse square of the 
distance. From (^61)- (51) we obtain 



R^ = 1.77 X lO^Aiy'AsAsmm^/^ cm . (29) 

At frequencies at which scattering is larger than absorption, the emergent 
spectrum from the last scattering surface is (Zeldovich and Shakura 1969) 
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where kjf = 3.69 x 10^T~^^'^i'~^g(l — e~^'^/^'^^ cin^ is the bremsstrahlung 
absorption coefficient, B^{T) is the Planck function, g is the Gaunt factor, T 
is the envelope temperature, and z/ is the frequency of the emergent emission. 
The shape of the spectrum determines the observed color temperature of the 
envelope. The temperature derived from Wien's displacement law (which is 
occasionally used instead of the color temperature) is mainly determined by 
the plasma temperature ; m and rh actually affect only the flux and shape 
while leaving the position of the peak in the spectrum unchanged (at those 
values at which the emergent spectrum is not a blackbody spectrum near the 
peak) (see ( pO| ) and Fig. ^). The frequency-integrated expression ( |30[ ) yields 
the flux Q = 5AR^ T^^/*erg/s cm^ (the Gaunt factor is set equal to 1) and 
is an analog of Stefan-Boltzmann's law (Sunyaev and Shakura 1974). The 
envelope luminosity is then 

L = AttRIq = 2.15 X 10^^ T^'^R]'^^ erg/s , (31) 

where T4 = T/IO"^ K and -R12 = -R^/10^^ cm. The departures from Stefan- 
Boltzmann's law become significant at T > 2.3 x 10'^-Ri2 K. Assuming that 
the spectrum has a nonblackbody shape (i.e., when scattering dominates), if 
ky < Aa"^ R , we find that this condition for the emission at A = 5500 A takes 
the form m > 186 tti"^/^ (see also Fig. |^). Using (^), we thus obtain 



L = 1.8 X 10^6(^1 JA2A3)'/^m^/^m2V8Tf/' erg/s . (32) 



SS 433 

Object no. 433 from the catalog of Stephenson and Sanduleak (1977) is fa- 
mous primarily for the two oppositely directed jets which move at a velocity 
of 0.26c (Fabian and Rees (1979). According to the current model, the jets 
are generated at the center of an accretion disk around one of the components 
of the close binary system, in which supercritical accretion takes place. Mass 
flow onto the compact object proceeds in the binary on the thermal time 
scale of a normal star at a rate of lO~^-lO~'^M0/yr (Cherepashchuk 1981; 
Van den Heuvel 1981). The mass loss rate in the envelope is 10~^-10~^MQ/yr 
(Shklovskii 1981; Dopita and Cherepashchuk 1981; Van den Heuvel 1981). 
Assuming a lower limit on the flow rate, we find from (||) that the following 
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relation holds for SS 433: 

mm ^ 4.7 X 10^^ , (33) 

where m = Mo/Mcr, and -^ is the accretion efficiency in a thin disk (bears 
no direct relation to the actual accretion efficiency, which is Ltot/iMoC^) = 

l/m{'dLtot/LEdd) )■ 

For our estimates we consider two cases: (1) ?7i = 10, m = 1000; and 
(2) m = 1 -^ 2, m = 10000. However, a black hole and a neutron star as the 
compact object represent markedly different situations, because, in contrast 
to a neutron star, a black hole has an event horizon under which the matter 
and energy accreted in an advective disk can freely fall. From (^) we derive 
the total luminosities of outflowing advective accretion disks: (1) Ltot ~ 
5.4 LEdd,BH ^ 7 X 10^^ erg/s (black hole) and (2) Ltot ~ 7.0 LEdd,NS ~ 
9 X 10^^ erg/s (neutron star). For supercritical accretion onto a neutron star, 
the bulk of the energy from the magnetospheric region can be carried away 
by neutrinos (Basko and Sunyaev 1976). 

Using ( p3D and (|29|) , we can find a relation between the uncertain mass 
of the compact object in SS 433 and the radius of the visible envelope (see 
Fig. I): 

Ri2 ^ 14.6 #/=^m-^/2 _ (34) 

Various theoretical and observational estimates yield an expanding-envelope 
radius of ~ 10^^-10^^ cm for SS 433 (see, e.g., Fabrika 1984; Bochkarev and 
Karitskaya 1985). 

In the envelope of SS 433 for rri = 10 and rh = 1000 or even for m = 1^2 
and rh = 10000 , electron scattering dominates over absorption (Fig. ^. The 
color temperature determined from the peak flux from SS 433 (Murdin et 
al. 1980) is ~ 25000 K, which corresponds to the plasma temperature in the 
envelope T = 3.5 x 10^ K (Fig. |p . From ( |52| ) we then derive L ^ 3 x 10^'' erg/s 
for m = 10 and L ^ 1 x 10^^ erg/s for m, = 2. This luminosity estimate 
depends on the mass flow rate in the envelope. The formula Voo{Rs) = 
0.62 c m~^/^ yields ~ 6000 km/s for IOMq. The observed Doppler widths of 
the "stationary" Ha and He II lines are ~ (1 — 2) x 10^ km/s. 
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Conclusion 

Here, we considered supercritical accretion onto a compact object. As a 
result, we obtained an analytic solution of the accretion- disk structure using 
classical equations. We developed a numerical scheme for calculating the 
accretion-disk structure with advection, i.e., with radial heat transport by 
the accreted matter. We obtained an estimate of the visible radius and 
luminosity of the envelope composed of outflowing matter as a function of 
the rate of accretion onto a black hole by taking into account bremsstrahlung 
absorption and Thomson scattering of photons. 

The physical quantities in the equations we considered were averaged 
over the height. Finding an exact solution is a difficult and interesting 
problem, which requires considering two-dimensional hydrodynamic equa- 
tions and which may also be associated with the problems of the generation 
of jets from sources with black holes. 

Here, we also assumed that the orbital particle velocity was described 
by Kepler's law. In this case, we ignored the contribution of pressure to 
the equation of motion (we took into account the pressure only when we 
calculated the vertical disk structure, Eq. (|)). In distant regions, at r >> 
Ro, this is quite acceptable. In Fig. 0, the speed of sound, the radial velocity 
over the disk, and the orbital Keplerian velocity JGM/r are plotted against 
radius. Near the sonic point, i.e., at the distance at which the speed of sound 
is equal to the radial velocity, a = w^, the contribution of pressure is large, 
and, as we see, the orbital and radial velocities are of the same order of 
magnitude there, i.e., the particle orbits cease to be Keplerian. 

The energy-wind equation differs from the relation that is used in the 
standard stellar-wind theory, 

TUpV ~ — . (35) 

If this law is used, then the efficiency of mass outflow from the disk decreases 
by several order of magnitude. Figure ^ shows the results of our numerical 
calculations for an advective disk with mass outflow that follows the law (|35|). 
It turns out that using this law generally reduces the outflow efficiency by 
more than two orders of magnitude and apparently cannot explain the rates 
of mass outflow from SS 433, which are comparable to the accretion rate in 
this system. 
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A distinctive feature of advective disks is that only a small fraction of 
the liberated gravitational energy of the accreted matter is emitted, while 
the heat is mainly transported in the disk. If the central star is a black hole, 
then matter together with energy fall into it without returning. By now the 
accretion-disk theory has progressed appreciably and has been applied to ob- 
served sources. It turns out that advection can be of considerable importance 
in supercritical accretion and at low accretion rates (see, e.g., Abramowicz 
et al. 1988, 1991; Narayan and Yi 1994; Narayan et al. 1997). In particu- 
lar, in the opinion of these authors, the model of a subcritical, optically thin, 
two-temperature advection-dominated disk around a supermassive black hole 
can explain the observations of low-luminosity galactic nuclei (Narayan et al. 
1995; Lasotaet al. 1996). 

Observations of SS 433 show that the disk loses up to 90% of the matter 
flowing to its outer boundary, producing an outflowing envelope around the 
compact star. Thus, the problem of calculating the structure of a noncon- 
servative accretion disk arises. In this paper, we considered an approximate 
solution of this problem. We found numerical solutions to this simplified 
problem (one-dimensional, etc.) with advection for a wide range of accretion 
rates above the critical value. We thus show that, in principle, substantial 
loss of the accreted matter from the accretion disk is possible, and calculate 
its value. 
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Figure 1: Relative thickness, accretion rate, fractions of emitted heat and 
advection versus radius. The sohd hne, the hne with short dashes, and the 
hne with long dashes represent an analytic solution with outflow, and advec- 
tive conservative disk, and an advective nonconservative disk, respectively. 
The distance is normalized to the inner disk radius; m = 1000. 
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Figure 2: Ratio of the time of photon diffusion to the surface to the charac- 
teristic time of radial matter motion in the disk versus radius. The numer- 
ical calculations were performed for a nonconservative disk with advection; 
rh =10, 100, 1000 and a = 0.5. 
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Figure 3: Change of the accretion rate in the disk due to mass outflow for 
various initial accretion rates. The curves (from left to right) represent the 
numerical calculations for an advective disk with m = 100, 500, 10'^, 10^. 
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Figure 4: Spectrum of the envelope: a blackbody spectrum at low frequen- 
cies; the emergent flux decreases with increasing frequency, when electron 
scattering begins to dominate over bremsstrahlung absorption. The spectra 
are given for m = 2, m = 5000 and m = 10, ih = 1000 (two upper close 
curves); m = 10®, m = 10 (middle curve); and m = 10®, m = 10^(lower 
curve). The plasma temperature in the envelope is T = 3.5 x 10'^ K. The 
temperature derived from Wien's displacement law is ~ 2.5 x lO'^K. 
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Figure 5: Curves that represent the lower boundaries of the domains of 
{m,m) values in which the optical depth for scattering is larger than the 
optical depth for absorption aX u = 10^^ Hz (lower curve), z/ = 4.5 x 10^^ Hz 
(middle curve), and z/ = 1.5 x 10^^ Hz (upper curve). 
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Figure 6: Visible envelope radius versus unknown mass of the compact 
component in SS 433. 
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Figure 7: Speed of sound a (solid line), radial velocity in the disk Vr (long 
dashes), and Keplerian velocity f<^ (dot-dashed line) versus radius. The nu- 
merical calculations were performed for a nonconservative disk with advection 
for m = 10, m = 1000, and a = 0.5. 
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Figure 8: Ratio of the accretion rate in the disk to the initial accretion 
rate versus radius for mass outflow that follows the law (0). The numerical 
calculations were performed for an advective disk with m = 10 and m = 
1000. 



26 



